Numerical study of a model for non-equilibrium wetting 
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We revisit the scaling properties of a model for non-equilibrium wetting [Phys. Rev. Lett. 79, 
2710 (1997)], correcting previous estimates of the critical exponents and providing a complete scaling 
scheme. Moreover, we investigate a special point in the phase diagram, where the model exhibits 
a roughening transition related to directed percolation. We argue that in the vicinity of this point 
evaporation from the middle of plateaus can be interpreted as an external field in the language of 
directed percolation. This analogy allows us to compute the crossover exponent and to predict the 
form of the phase transition line close to its terminal point. 
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I. INTRODUCTION 

While wetting of surfaces at or near thermal equilib- 
rium is well understood f]|, the study of wetting phe- 
nomena far from equilibrium is a challenging new field. 
In the past decade there have been numerous theoretical 
studies addressing the question whether non-equihbrium 
conditions may lead to different physical phenomena near 
the wetting transition. Most of these studies are based 
on particular lattice models [3, y, [j] or phenomenolog- 
ical Langevin equations [5|. They all have in common 
that the wetting layer is modelled by a non-equilibrium 
growth process of a d-dimensional interface combined 
with a hard-core wall which represents the surface of the 
substrate. 

The theoretical interest in non-equilibrium wetting 
stems from the fact that various scale-invariant prop- 
erties are found to be universal, i.e., they are dictated 
by the symmetries of the model irrespective of micro- 
scopic details. For a growth process without a substrate 
(free interface), the most prominent universality classes 
are the Edwards Wilkinson (EW) and the Kardar-Parisi- 
Zhang (KPZ) universality classes [y, Q- These classes 
describe the asymptotic scaling behavior of the rough- 
ening interface and are characterized by a certain set of 
exponents and scaling functions [q . In the corresponding 
Langevin equations, the KPZ class differs from the EW 
class by a non-linear term that breaks reflection symme- 
try in height direction. In experimental setups, where 
this symmetry is generally broken, KPZ behavior is ex- 
pected to be generic while linear growth (EW behavior) 
can be considered as a special case. 

Non-equilibrium wetting is usually modelled as a 
stochastic growth process on top of a hard-core substrate 
at height zero. Varying the growth rate the presence of a 
substrate induces a wetting transition from a bound to a 
moving phase. Concerning the scaling properties of the 
interface the substrate plays the role of a boundary: It 
does not change the universality class of the growth pro- 
cess itself, instead it imposes additional features. More 
specifically, it gives rise to an additional order parame- 



ter and an associated critical exponent. The situation 
is similar as e.g. in the Ising model, where a boundary 
induces an additional surface critical exponent. We will 
refer to these extended universality classes as the bounded 
Edwards- Wilkinson (bEW) and the bounded Kardar- 
Parisi-Zhang (bKPZ) class. For KPZ-type growth, how- 
ever, where the reflection symmetry is broken, it turns 
out that the new exponent also depends on the sign of 
the non-linear term. Therefore, one has to distinguish 
two different bounded KPZ classes, which we shall de- 
note by the acronyms bKPZ+ and bKPZ- according to 
the sign of the nonlinear term. 

One of these models for non-equilibrium wetting, which 
has been studied intensively in the past, was introduced 
a decade ago in Ref. ^. It is a restricted solid-on-solid 
(RSOS) growth process on a 1-dimensional lattice (see 
Sect. 2 for details), where the substrate is introduced by 
imposing the condition that all heights have to be non- 
negative. The phase diagram of this model is shown in 
Fig. [D Depending on the growth rate q and the evapora- 
tion rate p, the model exhibits a wetting transition from 
a bound to a moving phase. For p — 1 the transition was 
shown to belong to the bEW universality class, while for 
p ^ 1 the transition belongs to one of the two bounded 
KPZ classes. However, the reported estimates for the 
critical exponents are still contradictory. One aim of the 
present paper is to clarify this issue and to confirm KPZ 
scaling along the whole line except for p = and p — 1. 

As shown in the figure, the transition line ends at the 
left terminal point at p = and q^^ = 0.3993(1). In 
this point the wetting model reduces to a special growth 
process which was studied earlier by Alon et al. ^] and 
has the special property that evaporation from completed 
layers is forbidden. This means that the interface can- 
not have a negative growth velocity and therefore the 
presence of a wall makes no difference. Varying q while 
keeping p — Q the model displays a roughening tran- 
sition at g = <1?^ ■ It was shown that the dynam- 
ics of sites at the bottom layer can be related to a di- 
rected percolation (DP) [1^ process, which is another 
class of non-equilibrium phase transitions different from 
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FIG. 1: Phase diagram of the wetting model with r = 1 



both bEW and bKPZ. Extending this analogy, it was ar- 
gued that the dynamics of the first few layers may be de- 
scribed in terms of unidirectionally coupled growth pro- 
cesses [ii!,|i3- 

An open question, which will be addressed in the 
present work, concerns the crossover from DP to the 
bKPZ- class in the vicinity of the DP point. In order 
to describe this crossover, we propose to interpret evap- 
oration from the middle of a plateau with a small rate 
p <^ 1 as a weak external field in the language of DP. 
This allows us to express the crossover exponent, which 
determines the characteristic shape of the transition line 
as it approaches the DP point, in terms of the response 
exponent of directed percolation. 

The paper is organized as follows. In Section 2 we es- 
timate the critical exponents for the model for a wetting 
transition by off-critical, time-dependent, and finite-size 
simulations. Moreover, we summarize a scaling picture 
for the bKPZib and bEW classes. Section 3 is devoted 
to the crossover from DP class [p — 0) to bKPZ- class 
(0 < p < 1) in the vicinity of the left terminal point of the 
transition line. Finally we demonstrate that the proposed 
interpretation of the wetting model as a DP process with 
a external field is in agreement with numerical results. 



II. ESTIMATION OF THE CRITICAL 
EXPONENTS 

A. Definition of the model 

The model for non-equilibrium wetting proposed in 
Ref. [3| is defined on a one-dimensional lattice with L sites 
and periodic boundary conditions. Each site i is associ- 
ated with a variable hi — 0, 1, 2, 3... which describes the 
height of the interface at site i. The interface obeys the 



FIG. 2: Transition rates for the wetting model 

restricted solid-on-solid (RSOS) condition 

\hi - hi±i\ < 1, 



(1) 



i.e., the heights at neighboring sites may differ by at most 
one unit. 

The interface evolves in time by random-sequential up- 
dates as follows. For each update a site i of the lattice 
is randomly chosen and one of the following processes is 
selected (cf. Fig. O: 

(a) deposition of a particle {hi —> hi + I) with rate q, 

(b) evaporation of a particle (hi — > /ij — 1) at the edges 
of plateaus with rate r, 

(c) evaporation of a particle (hi -^ hi — I) from the 
middle of a plateau with rate p. 

A move is rejected if it would violate the RSOS condi- 
tion. Moreover, the substrate is introduced by imposing 
the restriction that evaporation at zero height is forbid- 
den. Each run starts with a flat interface at zero height. 
Without loss of generality we set r — I. 

The phase diagram of the model is shown in Fig. [T] 
The transition is controlled by the growth rate q. Above 
the critical line the model is in the wet or moving phase, 
where the interface roughens and propagates at constant 
velocity. Below the critical line the interface remains 
bound and fluctuates close to the wall. The value of 
the evaporation rate p determines the type of the phase 
transition. From the physical point of view small values 
of p are more realistic since evaporation at the edges of 
a plateau is usually more likely than in the middle. 

The order parameter for the wetting transition is the 
density of sites at zero height uq. In the bound phase 
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FIG. 3; Off-critical simulations. Density of sites with height 
zero no (left) and the interface width w (right) as functions 
of the distance from the critical point qc — 0.4295(1). With 
p — 0.001, L — 4096 and 100 independent realizations. 



near the transition line no goes to zero as 
no ~ {qc - qY , 



(2) 



where qdp) is the transition point and /3 is a critical 
exponent. Similarly, the interface width w diverges at 
the transition as 



w 



{qc - q) 



< 



(3) 



where C is another critical exponent. 

The vertical line p = is special in so far as a layer, 
once completed, cannot evaporate again. As already 
mentioned in the introduction, this special case was stud- 
ied in Ref. [9| and the exponent /3 was found to be in 
agreement with the DP exponent /3 = 0.27649(4) p^ . 
Moreover, it was shown numerically that for p = the 
interface width diverges logarithmically near the critical 
point jl4j | , which is consistent with the critical exponent 

C = o. 

Another special case is p = 1 , where the wetting tran- 
sition belongs to the bounded Edwards- Wilkinson (EW) 
class. In this case the stationary state in the bound phase 
can be computed exactly by transfer matrix methods [3] , 
which allows one to calculate the critical exponents /3 = 1 
and C = 1/3. 

For < p < 1 the model shows a different critical 
behavior as in the equilibrium case p = 1. The first es- 
timates reported in y, using off-critical numerical sim- 
ulations at p = 0.05, are /3 = 1.51(6) and C = 0.41(3). 
The purpose of this section is to revise these values, con- 
firming the conjecture that in this regime the transition 
belongs to the bKPZ- class. Similarly, for p > 1 the 
transition is expected to belong to the bKPZ+ class. 



B. Off-critical simulations 

First we calculate the exponents /3 and C by means of 
off-critical simulations. To reasons to be explained below, 



we use a very small value p 
shown in Fig. [3] we obtain 



0.001. From the graphs 



^=1.67(5), C = 0.41(5). 



(4) 



We believe that this estimate of (3 is larger than the one 
obtained previously in [4] because of a crossover from EW 
to KPZ behavior. This crossover is known to be notori- 
ously slow and may cause the impression as if the critical 
exponents depended continuously on p, varying from the 
EW exponent (3=1 and some values larger than one. 
However, it seems that such an estimate is just an effec- 
tive exponent measured in the crossover regime. As one 
approaches the critical line for p < 1 and increases the 
numerical effort the effective exponent grows and slowly 
converges to the /true' KPZ exponent. This crossover 
from EW to KPZ is expected to become more pronounced 
if we move away from the equilibrium case p = 1 . For in- 
stance, with simulations at p = 0.9 we would obtain the 
effective exponent j3 = 1.02(5). That is why we chose 
such a small value for p. 



C. Finite-size simulations 

According to the standard scaling picture of non- 
equilibrium phase transitions, the spatial correlation 
length ^± near the critical point diverges as 



^_L ^ {qc - qY 



(5) 



For p = 0, where the model exhibits DP behavior, one 
obtains the DP values v^^ « 1.10. For < p < 1, 
where the model is in the bKPZ- class, we observe that 
the critical point, where the velocity of the free inter- 
face is zero, varies strongly with the system size. For 
instance, at p = 0.001 for L = 128,4096 we found 
qc{L) = 0.425(1), 0.4295(1). Therefore, it is near at hand 
to postulate the relation 



qc{^) - qc{L) ^ L-^'"^ , 



(6) 



where the qc{oo) is the extrapolated value of the critical 
threshold. With L = 64, 128, 256, 512, 1024 we obtain 
qc{oo) = 0.4295(3) and i/_l = 1.00(3). The data leading 
to this results, with A = qc{oo) — qc{L) as a function 
of L, is shown in Fig. |4l We note that the extrapolated 
value qc{oo) = 0.4295(3) already coincides with the value 
gc(4096) = 0.4295(1) within error bars. 

For finite growing interfaces, after an initial transient, 
the correlation length ^j_ becomes of the same order as 
the system size. When this happens the interface width 
saturates. The value at which the interface width satu- 



rates depends on L and scales as [15 1 



W,at{L)^L'', 



(7) 



where Wsat denotes the saturation value of the interface 
width and a is the so-called roughness exponent. This 







FIG. 4: Finite-size simulations. The figure shows the differ- 
ence A between the finite-size criticai point qc{L) and the 
extrapoiated criticai point qdoo) = 0.4295(3) as a function of 
L for p = 0.001. 



FIG. 5: Time-dependent simulations. The order parameter 
no at the critical point Qc = 0.4295(1) as a function of the 
number of Monte Carlo steps t with p — 0.001, L = 8192 and 
800 independent realizations. 



relation is valid above the critical point and at the critical 
point for p ^ 0. In one dimension one has a — 1/2 for 
both the KPZ and the EW class 15]. Eqs. ©, ^, and 
(|3l) imply the relation 



C = ^i 



(8) 



From h'^ — 1.00(3) we obtain ( — 0.50(1), which is signif- 
icantly larger than the numerical estimate in the previous 
subsection. 



for the KPZ class is 2 = 3/2 15]. With z = 3/2 and our 
previous estimate of :/^ we obtain :/|| = 1.50(5). Together 
with e = 1.15(3) this means that /3 = 1.73(6). 

We note that both values f3 = 1.73(6) and C = 0.50(2), 
obtained from finite-size and time-dependent simula- 
tions, are larger than the results obtained from off-critical 
simulations /3 — 1.67(5) and ( — 0.41(5). We beheve that 
these discrepancies can be traced back to the fact that 
the EW-KPZ crossover is more influential in off-critical 
simulations. 



D. Time-dependent simulations 

Likewise, the temporal correlation length ^|| diverges 
close to criticality as 



?|| ^ {Qc - q)' 



(9) 



where vh is the temporal critical exponent. From this 
relation and Eq. ([2]) we can conclude that the bottom 
layer density tiq decays at criticality according to a power 
law 



with 



no r^t 



/3/^ll 



(10) 



(11) 



The exponent 9 was measured in Ref.|4|, simulating the 
so-called single step model which is known to belong to 
the bKPZ- class. Our result 9 = 1.15(3), which is shown 
in Fig. \5\ is compatible with their estimate 9 = 1.184(10). 

The dynamical exponent z = vn/iy^, defined by 



^w-i 



_L' 



(12) 



E. Scaling picture for equilibrium and 
non-equilibrium wetting 



Let us now summarize the scaling picture and the val- 
ues of the critical exponents. One-dimensional wetting 
models, as the one defined in Ref. [2|], are characterized by 
four independent critical exponents, namely, three bulk 
exponents a,i'±,i>^^ which take simple fractional values, 
and one exponent /3 associated with the order parameter 
no J whose value is only known numerically. These ex- 
ponents, together with related exponents C = ^^±0; and 
9 — /3/t^ii , are listed in Table [D The scaling hypothesis 
states that any quantities are invariant under the scaling 
transformation 



r->br, i -> b^'t, 



w —I- b°'w, 



(13) 



where 6 is a scaling factor and A = g — gc denotes the 
distance from criticality. 

The three bulk exponents can be determined as follows. 
For p > two of them, namely, the exponents a = C^/i'jl 
and z — v^jvn^ are just the well-known bulk exponents 



case 


a 


z 


J^± 


V\\ 


C 


e 


/3 


DP 





1.58 


1.10 


1.73 





0.159 


0.276 


bKPZ- 


1/2 


3/2 


1 


3/2 


1/2 


1.184(10) 


1.776(15) 


bEW 


1/2 


2 


2/3 


4/3 


1/3 


3/4 


1 


bKPZ+ 


1/2 


3/2 


1 


3/2 


1/2 


0.228(5) 


0.342(8) 



TABLE I: List of the critical exponents. Most of DP expo- 
nents in the frist line come from [I^l , the exceptions are a and 
C, that come form [1J|. The bEW and bKPZ exponents come 
from an exact result [3] and the best numerical estimatives 
for the exponent Q obtained in |j| respectively combined with 
the scaling picture developed here. 



of the EW or KPZ universality classes. To obtain the 
third exponent, let us consider the propagation velocity 
of a free interface. For p ^ 0, the interface velocity varies 
linearly with the distance from the critical line, i.e. 

vr^q^q^. (14) 

Since the velocity is the temporal derivative of the mean 
height, which has the same scaling dimension as the 
width ui, we expect that w ^ {q ~ ^c)?!! • With equa- 
tions (|12p and ([7]) we arrive at the scaling relation 



6 



(9-gc)-^/^^-") 



hence 



(15) 



(16) 



This argument is hand-waving as it uses the properties 
of the free interface to predict the scaling properties of 
the wetting transition, but it is in agreement with all 
numerical observations. For example, for the KPZ case 
we have z/j^ = 1, in agreement with our numerical result, 
and for the EW case one obtains v\^ = 2/3, in agreement 
with an exact calculation [J]. 

We note that the determination of f3 (or likewise 9) 
in the KPZ regime remains a numerically challenging 
task. Currently the most precise estimates come from 
the single-step model investigated in Q, reporting the 
values e = 1.184(10) for the bKPZ- and 9 = 0.228(5) 
for the bKPZ-|- classes. The estimates obtained in the 
present simulations are consistent but not as precise, in- 
dicating that KPZ behavior of the model introduced in 
Ref. [2| is not as 'clean' as in the single-step model. 



III. THE LIMIT p^O: INTERPRETATION AS A 
DP PROCESS IN AN EXTERNAL FIELD 



As mentioned in the preceding section, the non- 
equilibrium wetting model introduced in Ref. [4] includes 
a special casep = 0, where it exhibits a transition belong- 
ing to the directed percolation (DP) universality class. 
The crossover from bKPZ- to DP has not been studied 
so far. 



The case p = is special in so far that a layer, once 
completed, cannot evaporate again. This means that the 
process does not feel the hard-core wall any more; it may 
be removed without changing the properties of the model. 
It was shown in Ref ^] that the sites at the actual bot- 
tom layer may be interpreted as the active sites of a DP 
process. This mapping is exact without RSOS condition 
but it remains effectively valid when the RSOS condition 
is imposed. 

The key observation of the present work is that the pro- 
cess controlled by the parameter p, namely, evaporation 
from the middle of a plateau, corresponds to spontaneous 
creation of active sites in the language of DP. In DP such 
a spontaneous creation of activity is interpreted as an 
external field conjugate to the order parameter and the 
corresponding scaling laws are well understood. In this 
section we use this analogy to predict the properties of 
the crossover from DP to bKPZ- in the non-equilibrium 
wetting process. 



Mapping the contact process in an external 
field to non-equilibrium wetting 



To understand the mapping between non-equilibrium 
wetting and DP in more detail, let us first consider the 
contact process (CP) [l3[ in an external field h, which is 
defined by the following dynamical rules with random- 
sequential updates: 



01(10) 



1 — >0 


with rate 


1, 


— >1 


with rate 


h, 


1(11) 


with rate 


A/2. 



(17) 



For h ^ Q one retrieves the usual CP which exhibits a DP 
transition. For h > this transition is destroyed because 
the model does no longer have an absorbing state. 

Let us now compare this process with the unrestricted 
variant of the growth model introduced in [9|, which 
evolves according to the following dynamical rules: 





h^ >h, + l 


with rate 


1, 




h, — .0 


with rate 


h, 


h.,- 


-> min{hi, ft-^+i} 


with rate 


A/2, 


hi — 


-^ min{hi,hijf-i} 


with rate 


A/2 



(18) 



It is straightforward to verify that the variable rji = Sh^.o 
follows exactly the dynamical rules given in ([TT]) and that 
the order parameter of the transition is the density of 
sites at the bottom layer. 

For the restricted variant of the model introduced in 
the above mapping is no longer exact. However, we argue 
that an external field h can be introduced in the restricted 
model in an effective way by modifying the dynamical 
rules according to Fig. [Sj This modification is motivated 
as follows. In the unrestricted version the likelihood for 
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FIG. 7: The difference qc — q!?^ as a function of p. 



FIG. 6: Transition rates for the restricted model 



evaporation from the middle of a plateau does not de- 
pend on the actual configuration of the interface in the 
vicinity. To establish a similar independence in the re- 
stricted version, we introduce an additional evaporation 
rate h which is the same for all interface configurations 
that respect the RSOS condition after an evaporation 
event. As usual, evaporation at zero height is forbidden. 

Obviously, the dynamical rules listed in Fig. [6] can be 
related to the wetting model by identifying the paranre- 
ters 



1, 



X + h, p — h. 



(19) 



Since h > 0, this relation is valid only in the region 
p < r. However, the mapping between the two mod- 
els is not one-to-one because their transition rates are 
slightly different. While in the dynamical rules of Fig. [6] 
evaporation with one neighbor at the same height hap- 
pens with rate A/2 + h, the corresponding event in the 
wetting model takes place with rate X + h. We expect 
that this minor difference does not change the critical be- 
havior. Therefore, we conclude that the interpretation of 
evaporation in the middle of plateaus as an external field 
in the language of DP is still valid even in the restricted 
variant. 



B. Prediction of the form of the critical Hne in the 
limit p ^ 0. 

We now demonstrate that the conjecture presented 
above can be used to predict the form of the phase 
transition line in the vicinity of the DP transition point 
{p,q) = (0,0.3993(1)). To this end we consider the sta- 
tionary values of the order parameter ng along the line 



Using the parameter h which is related to q and p by 
q = l/[(gf ^)-i + h] and p = h/[{q^P)-^ + h] we define 
the field exponent 5 by the asymptotic power law 



no 



(21) 



which is expected to be valid for small h. With numerical 
simulations we obtain the value 5 — 0.121(5). This value 
differs from the DP exponent (3 = 0.27649(4) measured 
in vertical direction along the line p = 0. Therefore, ap- 
proaching the DP point from different directions we find 
two different exponents. The smaller one (with a slower 
decay of no) is expected to dominate all other directions 
except the vertical one. This allows us to conclude that 
in horizontal direction, i.e., for q = q^ and small values 
of p, the order parameter vanishes as 



n^r^p 



This equation together with equation Q implies 



^DP 



9c 



' p" 



(22) 



(23) 



q 



-.DP 



(1-p)- 



(20) 



where y = 5/ fi = 0.44(2) is a crossover exponent. 

The crossover exponent describes how the critical line 
approaches the DP point in the sense that it determines 
the critcal line concavity near the DP point. In fact, 
plotting qc — q^^ versus p in a double-logarithmic plot 
(see Fig. [7]) one obtains a straight line with the slope 
0.43(2), which coincides with the value of y. This means 
that the interpretation of evaporation from the middle of 
a plateau as an external field in the language of DP yields 
the correct crossover exponent describing the curvature 
of the phase transition line. 



C. Critical behavior of the first few layers 

So far we considered only the density of sites at the 
bottom layer uq. Similarly one can study the density 
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FIG. 8: TUk as a function of h with Ac = 4.30195(5) ^ for the 
restricted model with the external field. In this case L = 4096 
and 1000 is the number of independent realizations. 



of sites mt. whose heights are less or equal than k. As 
shown in [111 Il2l | , these order parameters vary in vertical 
direction (p = 0) by power laws m^ ^ {<£^ — q)^'' with 
individual exponents /?o,/?i,/32, . . ., where /?o = /5 is the 
ordinary density exponent of DP. In the same way we can 
now define the exponents 6k by 



ink 



P 



(24) 



where Sq is just the exponent 6 of the preceding subsec- 
tion. 

As an example, Fig. [8] shows numerical measurements 
of irik for the restricted model. The estimates of S^ are 
listed in Table |TT1 All these values are in fair agreement 
with the numerical estimates for the unidirectionally cou- 
pled DP reported in [1^ . 



IV. CONCLUSION 

In this paper we have presented improved estimates 
of the critical exponents for the non-equilibrium wetting 
model of Ref. Q in the parameter range < p < 1. We 
have derived relations between the exponents, in agree- 
ment with the numerical results, and used them to give a 
complete set of the exponents along the whole transition 
line (see Table U. 

This work is focused on the case < p < 1, where 
the wetting transition belongs to the bKPZ- class. For 



model 


So 


Si 


S2 


unrestricted 


0.107(2) 


0.040(3) 


0.014(2) 


restricted 


0.104(5) 


0.039(5) 


0.008(5) 


wetting 


0.121(5) 


0.047(5) 


0.009(5) 



TABLE II: Numerical estimates of the exponents 5^ 

p > 1, where the transition belongs to the bKPZ-|- class, 
we could do a similar analysis but we would have to use 
extremely large values of p to overcome crossover effects, 
but numerical simulations turn out to be inefficient in 
this limit. 

Presently the most reliable results for the critical expo- 
nent 9 of the bounded KPZ classes were measured using 
the so-called single step model fij. This model has a 
moving wall and is always exactly at the critical point. 
However, the single step model does not allow one to 
perform off-critical simulations. This is the reason why 
the wetting model is more suitable to confirm the scaling 
relations discussed in this paper. 

The non-equilibrium wetting model introduced in [3| 
has another interesting feature, namely, a special transi- 
tion point at p == 0, where the critical behavior belongs to 
the directed percolation universality class. Surprisingly, 
even for very small p, the critical behavior changes en- 
tirely. For example, the exponent /3, which describes the 
density of sites at zero height, jumps from the DP value 
P w 0.28 for p = to a large value /? w 1.78 of the bKPZ- 
class. To our knowledge, this is the only case where the 
order parameter exponent /3 is larger than one. This is 
due to the fact that the interface at criticality has only 
very few contact points where it touches the wall. This 
number is so small that the correlation length ^_l of the 
interface fluctuations becomes smaller than the average 
distance between two contact points [4|, which is equal 
to the inverse of the density of sites at the bottom layer. 

The main result of this work is the conjecture that a 
wetting process for small p can be interpreted as a DP 
process in an external field. This conjecture allows us to 
interpret the crossover from the DP class to the bKPZ- 
class as a small external field h that eliminates the DP 
transition. Moreover, it is consistent with the scaling 
picture and can be confirmed by numerical simulations. 
Calculating the corresponding crossover exponent y — 
6/(3, we can predict the shape of the critical line near 
p = 0. 
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